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ABSTRACT 

This work describes a new instrument optimized for a detection of the neutral hydrogen 21cm 
power spectrum between redshifts of 0.5 — 1.5: the Baryon Acoustic Oscihation Broadband and 
Broad-beam (BAOBAB) Array. BAOBAB wih build on the efforts of a first generation of 21cm 
experiments which are targeting a detection of the signal from the Epoch of Reionization at z ~ 10. 
At z ~ 1, the emission from neutral hydrogen in self-shielded overdense halos also presents an accessible 
signal, since the dominant, synchrotron foreground emission is considerably fainter than at redshift 
10. The principle science driver for these observations are Baryon Acoustic Oscillations in the matter 
power spectrum which have the potential to act as a standard ruler and constrain the nature of dark 
energy. BAOBAB will fully correlate dual-polarization antenna tiles over the 600-900MHz band with 
a frequency resolution of 300 kHz and a system temperature of 50K. The number of antennas will grow 
in staged deployments, and reconfigurations of the array will allow for both traditional imaging and 
high power spectrum sensitivity operations. We present calculations of the power spectrum sensitivity 
for various array sizes, with a 35-element array measuring the cosmic neutral hydrogen fraction as a 
function of redshift, and a 132-element system detecting the BAO features in the power spectrum, 
yielding a 1.8% error on the z ^ 1 distance scale, and, in turn, significant improvements to constraints 
on the dark energy equation of state over an unprecedented range of redshifts from ~ 0.5 — 1.5. 
Subject headings: instrumentation: interferometers — cosmological parameters — distance scale — 
techniques: interferometric — large-scale structure of the universe 



1. INTRODUCTION 



The Baryon Acoustic Oscillation (BAO) features in the 
large-scale matter distribution have recently drawn at- 
tention as a standard ruler by whi ch the geometry of the 
unive rse can be directly measured ( Eisenstein et al.|1998 
1999). These features in the cosmic microwave back- 



ground (CMB) power spectrum and the matter power 
spectrum today are imprints from the acoustic oscilla- 
tions in the primordial photon-baryon plasma that re- 
combined at z ~ 1100. The features in the power spec- 
trum appear at multiples of the sound horizon scale at 
recombination, making them effective standard rulers. 
Measuring the BAO wiggles at several redshifts yields ge- 
ometric measurements of the universe — the Hubble pa- 
rameter, H{z), and the angular diameter distance, dA{z) 
— that constrain properties of the dark energy that dom- 
inates the cosmic energy content at z = and is the 
current leading theory for the accelerated expansion of 
thejiniverse. Since the first detection of the BAO signal 
( [Eisenstein et al.||2005| , several experiments have been 
undertaken to use these features for precision cosmology, 
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including the Sloan Digital Sky Survey HI Ba ryon Os 
dilation Spectroscopic Sur vey (SDSS-III BOSS; iS c hlegel 
et al l[2 QQ9p] WiggleZ ( |Drinkwater et al.| |201oV ^| and 



the Hobby- T?berly Telescope Da rk Energy Experiment 



(HETDEX; [Adams et al]|2011p j as well as a num 



ber of planned future experiment s, such as the Su baru 

"ar|[2 QT2l), Eu- 

'Schle 



Prim e Focus Spectro^aph ( PES; lEllis et 



did dAmendola et al.||2QT?), BigBOSS ("Schlegd et al. 



2011), and the Wide-Field Infrared Survey Telescope 
(WFIRSTp^ All of these experiments target individual 
galaxies witn spectroscopic observations. 

Rather than targeting individual objects, a 21cm inten- 
sity mapping experiment can detect fluctuations in neu- 
tral hydrogen emission on large scales ( Chang et al.|2008 



iWyithe et alT2QM| [Morales fc Wyith e 2010; Prit chard &: 
[Loeb 2011 ), with two dimensions corresponding to angles 
on the sky, and the third line-of-sight dimension aris- 
ing from the differential redshifting of 21cm line emis- 
sion as a function of distance. After reionization, the 
power spectrum of 21cm fluctuations is expected to be 
a biased tracer of the matter power-spectrum, since the 
remaining neutral hydrogen resides in high-density, self- 
shield ed regions such as in ga l axies and other coll apsed 



halos (Barkana & Loeb 2007 Madau et al. 1997). As 



a result, 21cm intensity mapping experiments present a 
promising complement to spectroscopic galaxy surveys 
for BAO science. Several 21cm intensity mapping ex- 
periments have been proposed, including the prototype 

http:/ /www. sdss3.org/surveys/boss.php/ 
http:/ /wigglez. swin.edu.au/site/ 
^'^ http://hetdex.org/ 

http://wfirst.gsfc.nasa.gov/ 
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Cylindrical Radio Telescope (CRT; fo rmerly HSHS, iPe] 

the Canadian 



terson et al. 


2006 




Seo et al. 


2010|i4 



Hydrogen Intensity Mapping Experiment (CHIME 
BAORadio ( Ansari et al. 2012a^b ), the BAO from Inte- 
grated Neutral Ga s Observations experiment (BINGO; 
Battye et al.||20l'2l and an on going experirnent w ith the 
Green Bank Telescope (GBT; [Chang et al |[2QTq1 ). 

The flexibility in angular and spectral responses of 
radio interferometers, which measure the power spec- 
trum both parallel and perpendicular to the line of sight 
( [Morales 200 5), gives 21cm BAO the ability to survey 
larger cosmological volumes and operate over a wider 
range of redshifts than current spectroscopic galaxy red- 
shift surveys. As a result, the 21cm BAO signal has 
the potential for probing expansion throughout and be- 
yond the critical epoch when dark energy comes to dom- 
inate the energy density of the universe. Furthermore, 
a 21cm intensity mapping experiment can probe red- 
shifts z > 0.5 with roughly uniform sensitivity, with- 
out complications arising from sky emission lines in the 
optical/near-infrared. The 21cm signal can be used to 
constrain the location of the BAO peaks as a function of 
redshift, and thereby measure the magnitude and time- 
evolution of dark energy. 

21cm BAO experiments can draw on the considerable 
investments in low-frequency radio astronomy developed 
in the past decade for studies of the Epoch of Reioniza- 
tion (EoR). In this paper, we present the Baryon Acous- 
tic Oscillation Broadband and Broad-beam (BAOBAB) 
array, a new experiment, building on the legacy of the 
Precision Array for Probi ng the Epoch of Reionization 
(PAPER; 'Parsons et al.| |2010|F^ and the Murchison 
Widefield Array (MWA; [Lonsdale et al.^ 2009pl for mea- 
suring the 21cm HI power spectrum at a redshiTt of 2: ~ 1. 
This paper is structured as follows: in ^ we present a 
system architecture for the BAOBAB instrument. In ^ 
we forecast the sensitivity and cosmological constraints 
that will be achieved by BAOBAB. We consider several 
possible challenges and extensions for this approach in 
Q and conclude in ^ Throughout this work we as- 
sume the WMAP7 best fit ACDM cosmological model: 
h 0.7, fl M = 0.27, Jlh = 0.046, ^de = 0.73, and 
(78 0.8 (Lar son et a L| [2011 ). 



2. THE BAO BROADBAND AND BROAD-BEAM ARRAY 

The past decade has seen significant progress in the 
design, construction, and calibration of low-frequency 
interferometric arrays toward the goal of detecting the 
highly-redshifted 21cm signal from the Epoch of Reion- 
ization. The technologies used in BAOBAB inherit from 
two EoR experiments — PAPER and the MWA — but 
with several significant modifications to optimize the in- 
strument for BAO science. The entire signal chain will 
be re-tuned to operate between 600-900 MHz. These 
frequencies corresponds to redshifted 21cm emission be- 
tween z = 0.58 — 1.37, a band chosen for several rea- 
sons. First, these moderate redshifts complement the 
undergoing lower redshift galaxy surveys like BOSS by 
probing what is currently a relatively unexplored vol- 



■^'^ http : //cmb . physics . wise . edu/people/lewis/webpage/ 

http: / /www. mcgihcosmology.ca/ chime 

http://eor.berkeley.edu/ 

http://www.mwatelescope.org/ 
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Figure 1. System diagram of the BAOBAB interferometer. 
Dual-polarization antenna signals at -103dBm enter an uncooled 
low-noise amplifier (LNA) with +12.5dB gain and a noise figure of 
0.4dB (30K). Second-stage amplifiers add +36dB gain before trans- 
mission through 30 meters of LMR400 50^2 cable to a central enclo- 
sure. The signal is bandpass filtered (600-900MHz) and amplified 
+40dB to the optimal -22dBm input level for the ADCs. Each 
antenna signal is digitized and channelized in ROACH F-engines, 
reordered in transmission through a 10Gb Ethernet switch, and 
sent to GPU GTX580 X-engines for cross-correlation. Raw visi- 
bility data are passed to a host computer and written to a RAID 
storage unit in the MIRIAD UV file format for post-processing. 

ume of the universe. Secondly, at these frequencies, com- 
mercially available amplifiers and cables provide suitable 
low-noise performance that would not be obtainable at 
higher frequencies (lower redshifts). Furthermore, this 
band also avoids the bright sky noise and ionospheric ef- 
fects that complicate lower frequency (higher redshift) 
observations. 

BAOBAB will be a non-phase-tracking, broadband ar- 
ray of beam-formed tiles. By lowering the field-of-view 
of each element, the use of beamformed tiles like the 
MWA will significantly increase BAOBAB 's power spec- 
trum sensitivity over that of an array of the equivalent 
number of single dipoles, without increasing the correla- 
tor demands. Each tile will consist of 4 scaled versions 
of the PAPER sleeved-dipole design and groundscreens, 
electronically beamformed to point to zenith. Two linear 
polarization signals from each tile enter a digital signal 
processor that computes both auto- and cross-correlation 
products and outputs the results locally to disk. A block 
diagram is given in Fig. Tj the key properties of the 
BAOBAB system are listed in Table [l] and are described 
in more detail in the remainder of this section. 



Table 1 

Proposed BAOBAB Array 



Operating Bandwidth 


600-900 MHz 


Number of Tiles 


32-128 


Collecting Area per Element 


2.6 m^ 


Gain per Element 


18 dBi 


Field-of-View 


0.045 str 


Receiver Noise Temperature 


40 K 


System Temperature 


50 K 


Maximum Imaging Baseline 


60 m 


Redundant Baseline Scale 


1.6 m 


^minj ^max 


0.01, 2.5 /iMpc-i 


Array Configuration 


Reconfigurable (see Fig. ^ 


Frequency Resolution 


300 kHz 


Snapshot Integration Time 


10 s 



2.1. 



Siting 

Since radio-frequency interference (RFI) is prevalent at 
these frequencies, BAOBAB will need to be located at a 
radio quiet site. The bottom panel of Figure [2] shows pre- 
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Figure 2. Top: Leuschner Observatory, with a prototype 2- 
element BAOBAB interferometer deployed. This system was devel- 
oped and deployed by students as part of a Fundamentals of Radio 
Astronomy class at UC Berkeley. Bottom: solar fringes measured 
with the BAOBAB-2 prototype at Leuschner. 

liminary measurements made by a prototype 2-antenna 
BAOBAB interferometer deployed at the Leuschner Ob- 
servatory near Berkeley, CA. At this site, only 40 MHz of 
a 400-800 MHz operating band show solar fringes uncor- 
rupted by RFI, demonstrating the need for the primary 
BAOBAB deployment to be located at a quieter site, 
such as the NRAO site near Green Bank, WV. Next- 
generation activities may take place at the Square Kilo- 
meter Array South Africa (SKA-SA) reserve in the Karoo 
desert. This site is currently occupied by the PAPER and 
MeerKAT arrays, and has been shown to be a pristine 
RFI environment ( Jacobs et al.|[2QllJ . 



2.2. Analog System 

With the drastic reduction in sky noise relative to 
EoR frequencies, BAOBAB's system temperature will 
be dominated by the analog electronics. These com- 
ponents must therefore be optimized to reduce receiver 
noise while maintaining the smooth spatial and spectral 
responses that are a hallmark of the PAPER design and 
a key component of the del ay spectrum for eground isola- 
tion approach presented in ~ 



after PI 2b) and discussed in 



Parsons et al 



([2012b) (here- 
analog system 

will include the collecting element (consisting of 4 an- 
tennas and reflectors), low- noise amplifier, coaxial cable, 
and receiver. 

The BAOBA B element will begin with a 1/5-scale PA- 
PER antenna (Parsons et al. 2010), as shown in Figure 
|3l This desi gn is a dual-polarized version of the sleeved 
dipole design that uses a twin-resonance structure con- 
sisting of a pair of crossed dipoles located between a pair 
of thin aluminum disks. The element's reliability has 





Figure 3. Top: A prototype BAOBAB dipole antenna, designed 
as a 1/5 scale model of a PAPER dipole. Bottom: BAOBAB tile 
design with 4 dipoles and individual ground-screens. 

been demonstrated in PAPER arrays over the past sev- 
eral years. A trough refiector under each dipole will be 
used to increase the directivity toward zenith. The elec- 
tromagnetic behavior of the element was modeled ex- 
tensively for PAPER using GST Microwave Studio, and 
shown to perfor m as desired through calibration with ce- 



lestial sources in Pober et al. (2012|. The geometrically 
re-tuned prototype shown in the top panel of Figure |3] 
will be optimized to operate efficiently over the 600-900 
MHz band. 

Rather than deploy single elements like PAPER, 
BAOBAB will use a 2 x 2 tile of dipoles and ground- 
screens, as shown in Figure [3j A fixed zenith beam- 
former will be used to combine the signals, increasing 
the gain by 6 dB and reducing the field-of-view by a fac- 
tor of four. Both analog and digital beamformers are 
being investigated. A key issue is the mutual coupling, 
which should be reduced by the additional groundscreens 
between dipoles. The net effect is that for a fixed corre- 
lator size, the power- spec trum sensitivity is increased by 



a factor of four (see ^3.2). 



The amplifier designed for PAPER has a measured 



the 120-170 MHz band 
cation to BAO at z ~ 



(|Parsons et al. 


20101) 


1, we will moc 


lify tt 



For appli- 
lis amplifier 



design to operate from 600-900 MHz. Besides re-tuning 
the filter and amplifier circuits, however, one of the major 
activities in this modification will be to reduce the noise 
temperature of the front-end amplifier in order to obtain 
a target system temperature of 50 K. This change refiects 
one of the key differences between the BAO and EoR 
foregrounds. System noise in the EoR band is dominated 
by ~300 K sky noise from galactic synchrotron emission. 
In the BAO band, the sky temperature is reduced to ~10 
K, making the front-end amplifier the leading source of 
noise. Uncooled commercial UHF-band amplifier tran- 
sistors based on GasFET or HEMT technology can re- 
liably achieve noise figures of 0.4 dB, corresponding to 
a receiver temperature of 30K. A prototype BAOBAB 
balun/amplifier using a Hittite HMG617LP3 LNA with 
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Figure 4. The prototype balun and Hittite HMC617LP3 LNA for 
BAOBAB. The amphfier adds +30dB of gain with a quoted noise 
figure of 0.5 dB. 




Figure 5. PAPER'S 128-input correlator (shown) fohows the 
packetized frequency-cross-multiply (FX) architecture developed 
by the Center for Astronomy Signal Processing and Electronics 
Research (CASPER). Shown are 16 ROACH F-engines (left) and 
2 dual-CPU box X-engines (right). The first-generation BAOBAB 
correlator will modify the 64-input, 100-MHz PAPER correlator 
to become a 32-input, 300-MHz correlator. It will employ eight 
ROACH boards for spectral processing and four dual-CPU boxes 
for cross-multiplication. A 10-CbE switch is used to route data 
between boards. 

a quoted noise figure of 0.5 dB is shown in Figure [4j tests 
are underway to determine the noise temperature of the 
complete system. 

2.3. Digital System 

The BAOBAB correlator will follow the scalable corre- 
lator design used by PAPER and other members of the 
international Collaboration for Astronomy Signal Pro- 
cessing and Electronics Research (CASPERl^ a real- 
time digital correlator based on Field-Programmable 
Gate Array (F PGA) processors an d graphics processing 
units (CPUs) dParsons et "aI][2QQ8> |Clark et al. 2 011). 
The correlator architecture we employ uses modular sig- 
nal processing electronics and packetized communication 
protocols to build correlators that are flexible in the 
number of antennas correlated and the bandwidth cor- 
related by each antenna. A photograph of a 128-input 
FPGA/GPU correlator is shown in Fig. [s] 

The generic FX correlator architecture we implement 
consists of modules responsible for digitizing and chan- 

https://casper.berkeley.edu 



nelizing each set of antenna inputs (F-Engines), followed 
by a set of signal processing modules responsible for 
cross-multiplying all antennas and polarizations for a sin- 
gle frequency (X-Engines) and accumulating the results. 
Unique to this architecture, signal processing engines 
transmit packetized data through commercially available 
10 Gb Ethernet switches that are responsible for rout- 
ing data between boards. This architecture, along with 
analog-to-digital converters, modular FPGA-based sig- 
nal processing engines, and a software environment for 
programming, debugging, and running them, were devel- 
oped in collaboration with CASPER at the University 



of California, Berkeley (Parsons et al. 2008). The flexi- 
bility and modularity oi this correlator design shortens 
development time, and in this case, allows an existing 
64-input, 100-MHz PAPER correlator with 8 ROACH- 
boards and 4 dual- GPU boxes to be straight-forwardly 
modified to become a 32-input, 300-MHz BAOBAB cor- 
relator using the same boards and signal processing li- 
braries. A forthcoming publication on this correlator will 
be presented in Ali et al. (2013). 

2.4. Configuration 

BAOBAB will employ small antennas and above- 
ground cabling with relatively inexpensive LMR400 50- 
Ohm coaxial cables; these cables will not be buried, al- 
lowing BAOBAB to easily change between different array 
configurations by moving antenna elements. Following 
the principles outlined in (P arsons et al., 2012a, hereafter 
P12a), BAOBAB will employ a minimum-redundancy 
imaging configurations for characterizing foregrounds 
with minimal sidelobes and maximum-redundancy con- 
figurations to repeatedly sample the same locations in 
the uv-plane, substantially improving sensitivity to the 
three-dimensional power spectrum of 21cm emission at 
z ^ 1. Although future experiments may target a range 
of angular scales to map 21cm emission in the plane of the 
sky, by focusing on a limited number of Fourier modes, 
these maximum-redundancy configurations can improve 
sensitivity to the power spectrum by an order of magni- 
tude or more in mK^, relative to an equivalent observa- 
tion with a minimum-redundancy configurati on. 

However, as will be discussed further in ^ 3.2.1[ the 
mapping of baseline length to a transverse /c-mode on the 
sky is significantly larger for BAOBAB than for PAPER. 
In order to probe the relatively large-scale BAO features, 
then, BAOBAB will use the most compact configurations 
possible for its power spectrum measurements. Such an 
array configuration for a 35-element system, as well as 
that of a 32-element imaging configuration, are shown 
in Figure [6) The tiles are spaced 1.6m apart, effectively 
touching end-to-end. Investigations of cross-talk and mu- 
tual coupling will take place during an early prototype 
of the system; it may be the case that a phase-switch or 
additional shielding between tiles will be necessary to ac- 
commodate the short baselines required by BAO science. 

With a modular CASPER correlator increasingly able 
to process larger numbers of antenna inputs, BAOBAB 
naturally lends itself to a staged approach. Early < 16 
tile-element prototypes will characterize system perfor- 
mance, while a subsequent ~ 32-element array will study 
foreground emission and constrain the neutral hydrogen 
fraction as a function of redshift with a measurement of 
the 21cm power spectrum (p.4[). A ~ 128-tile version of 
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Figure 6. BAOBAB array configurations, plotted in meters east- 
west (horizontal axis) and north-south (vertical axis). BAOBAB 
will use above-ground cabling to allow antennas to be moved 
into a minimum-redundancy configuration (left) for imaging fore- 
grounds, or a maximum-redundancy configuration (right) for en- 
hanced power-spectrum sensitivity. Each square represents one 
tile, and each "+" one dipole. 

BAOBAB will measure BAO features and provide sub- 
stantial improvements over our current constraints on 
t he eq uation of state and time evolution of dark energy 



3. PREDICTED COSMOLOGICAL CONSTRAINTS FROM 
BAOBAB 

In this section we present predictions for forthcom- 
ing cosmological constraints for several iterations of the 
BAOBAB instrument. We begin by reviewing the pre- 
dicted signal strength for the cosmological 21cm power 
spectrum in §3.1 J In ^3^, we adapt the power spectrum 
sensitivity calculations of PI 2a for an array operating 
at z ~ 1, including the effects of sample variance and 
shot noise. In |3.3| we briefly review the delay spectrum 
foreground removal procedure of PI 2b. While a detailed 
study of foregrounds is beyond the scope of this paper, it 
is worthwhile to discuss the implications of the technique 
on which Fourier modes of the 21cm power spectrum 
will be accessible. We conclude the section by presenting 
forecasts for a 35- and 132-element BAOBAB system in 
§3.4| and ^3^ respectively, including Fisher matrix pre- 
dictions for constraints on the dark energy equation of 
state in the latter section. In the discussion of Q we 
explore possible directions for improvement with larger 
BAOBAB arrays. 

3.1. The 21cm Power Spectrum 

As with galaxy redshift surveys, a 3D map of the neu- 
tral hydrogen in the universe will serve as a tracer of the 
underlying dark matter power spectrum. The brightness 
of the observable radio 21cm signal will depend on the 
cosmological neutral hydrogen fraction, as well as the 
bias of hydrogen cont aining haloes with respect to the 
dark matter (Barkana & Loeb||2007l iMadau et al. 111997 



Ansari et al.||2Q12b| ): 



(1) 



f2i{z) 0.084mK 



{l^z)^h 



fm{z) 



0.044 0.01 



(2) 



where T2i{z) is the mean 21cm brightness temperature 
at redshift z; P{k,z) is the linear matter power spec- 
trum; b is the bias factor of HI containing halos with 



a. 
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Figure 7. Top: The predicted real space (i.e., /i = 0) power 
spectrum of 21cm emission in our fiducial model at three redshifts: 
0.67, 0.89 and 1.18 (corresponding to frequencies of 850, 750 and 
650 MHz). The z = 0.89 and z = 1.18 power spectra overlap, and, 
effectively, the signal is at the same strength for all three redshifts. 
Matter power spectrum predictions come from CAMB. Bottom: 

The equivalent dimensionless power spectra, A'^(k) = ^^P(k). 

respect to the dark matter; fm{z) is the mass fraction of 
neutral hydrogen with respect to the overall cosmological 
baryon content (i.e., = f^b)] is the cosmological 
constant, and Qrn and Qb are the matter and baryon 
density in units of the critical density, respectively. 

We plot the predicted 21cm brightness temperature 
power spectrum, for our fiducial model at red- 

shifts of 0.67, 0.89 and 1.18 (corresponding to frequen- 
cies of 850, 750, and 650 MHz, respectively) in Figure 
[7[ Predi ctions f or the matter power spectrum come from 
CAMB ( [Lewis et al. 2000 We also plot the dimension- 
less power spectrum, A^(/c) = ^P{k). For the remain- 
der of this paper will we primarily express our results in 
terms of A^^^(/c), due to its more intuitive units (mK^) 
and simple physical interpretation as the variance per 
logarithmic k bin. 

In practice, the 21cm power spectrum is measured in 
redshift space, which can b e, at linear o rder, related to 
the real-space PT^iik^z) as ( [Kaiser [p. 98 7[ ): 

(3) 



PT,,(fc,^) = (l + /?M')'Pr.,(fc,2) 



where /a = k • z is the wavevector k projected along the 
line-of-sight i, and P = f{Q)/b from linear theory where 
f{rt) ^ rtjn{z)^ is the dimensionless linear growth rate 
and 7 = 0.557 for ACDM cosmologies. To incorporate 

http://camb.info/ 
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the effects of redshift-space distortions in the fohowing 
calculations, we reduce our thermal noise error bars in 
/c-space by the factor of (1 + however, we do not 

attempt to constrain cosmological parameters like /3 by 
measuring the power spectrum as a function of /i. We 
assume fiducial values of /hi = 0.015 and 6 = 0.75, chosen 
to agree with fmb = 0.012 ± 0.003 as measured by 
Chang et al. (2010). To obtain individual constraints on 
these parameters, one will need to measure redshift-space 
distortions themselves. 

3.2. Sensitivity of an Array to the 21cm Signal 

There are three independent sources of statistical un- 
certainty in a 21cm power spectrum measurement: ther- 
mal noise in the inter ferometric visibilities, sample vari- 
ance, and shot noise, of which the last is in some sense 
"signal," but still inhibits measurements of cosmologi- 
cal parameters. For the first-generation of 21cm experi- 
ments, thermal noise will be the dominant source of un- 
certainty in measurements of the power spectrum. We 
therefore calculate the effects of thermal noise indepen- 
de ntly in §3.2. 1[ We ad d in th e effects of sample variance 
in ^ 3 2.2| and argue in ^ 3.2.3 that shot noise can be ne- 
glected for these experiments. 

3.2.1. Thermal Noise 

Thermal noise, in addition to being the dominant 
source of uncertainty in first generation 21cm BAO ex- 
periments, is also likely to be least familiar to those used 
to optical redshift surveys. Given the limited collecting 
area of early experiments, reducing thermal noise con- 
tributions is of paramount importance, even at the ex- 
pense of the number of Fourier modes measured. Much 
of the work in this section builds off of the derivation 
of an interferometric array's power spectrum sensitivity 
presented in P12a. We begin with a version of their equa- 
tion 16, which gives the power spectrum of the thermal 
noise obtained from one integration of a dual-polarization 
baseline: 

^l^k)^X'Y^^^Tl,, (4) 

where X^Y is a scalar translating observed units to cos- 
mological distances in /i~^Mpc (X converts from angles 
on the sky to transverse distances, Y from bandwidth to 
line-of-sight distance), ^ is the solid angle of the primary 
beam of one element in steradians, and Tgys is the system 
temperature. 

The science goals of a BAO experiment are to actu- 
ally measure X and Y] that is, since the exact values of 
X and Y depend on the underlying cosmology, we can 
combine the known physical scale of BAO with the angu- 
lar and frequency scales in the observed signal to extract 
the detailed expansion history of the universe. For the 
purpose of a sensitivity derivation, however, the behav- 
ior of X and Y can be considered well enough known to 
compute fiducial values. 

X is related to the angular size distance, I^a, as 



X = D^{l^z)= I 
Jo 



c dz 



(5) 



with H{z) in the matter /dark-energy dom- 
inated epoch being approximately given by 



H{z) = i^o\/^M(l + + ^A- Numerical integra- 
tion for a flat universe with = 0.27, I^a = 0.73, and 
= 70 yields Da{z = 1) ^ 1680 proper Mpc (|Wright] 
2006). Ignoring the evolution in the angular diameter 
distance around 2: ~ 1, we can write: 



X ^ 1700(1 + z) 



Mpc 
radian * 



(6) 



Note that although we use this (admittedly somewhat 
poor) approximation for simplicity in deriving the rela- 
tions for thermal noise power spectra in this subsection, 
all the subsequent results include the full z-dependence 
of the angular diameter distance. 

A few more words are warranted concerning the mag- 
nitude of X at z = 1. Given the scaling of equation [5j 
a 16A-baseline (6.8m a z = 1) corresponds to a trans- 
verse wavemode of k± = 0.042 /iMpc~^, a non- negligible 
value compared to the first BAO peak at ^ 0.08 hMpc~^. 
Therefore, baselines longer than ~ 32A will lose access 
to the first peak and be less effective probes of cosmol- 
ogy, regardless of foreground effects to be discussed later. 
This scaling motivates the extremely compact configura- 
tions proposed for BAOBAB in ^ despite the possible 
systematics associated with such short baselinesp^ 

To compute the scaling between frequency, and co- 
moving line-of-sight distance, rios, we use 



dr] 



los 



cdz 
H{z)- 



(7) 



Since v{l z) = we have that dz/{l + z) = —dvjv 
so 

_ curios _ c(l + zf 



Y - - 

dv h'2iH{z) 

Evaluating the above numerically, we get 

(l + z)2 Mpc 



y = 3.0 



V^^m(1 + ^)^ + ^a MHz* 



(8) 



(9) 



Finally, we compute the product X^Y used in equation 



x^r ^ 2.93 



yn^^TT^pTOv sr • Hz 



(10) 



Nominally, X^F = 2S{h-^Mpc)^ st'^ Rz'^ at z = 10 
The other values in equation [4] are system-dependent 
parameters. The BAOBAB tiles have a considerably 
sized primary beam on the sky, ft ^ 150sq.deg ^ 0.045sr. 
However, this beam is significantly narrower than that of 
a single dipole, so that the use of beamforming results 

Epoch of Reionization experiments at z = 9, however, do not 
find themselves hmited by the transverse modes probed. X{z = 



9) 



9360 



Mpc 



(P12a), so that a 16A-basehne corresponds to 



radian 

k± = 0.015 hMpc~^ . With this scahng, the effect of the k± compo- 
nent on the measured power spectrum wih always be sub-dominant 
to the foreground effects described in PI 2b. 

The magnitude of X^Y at 2; = 1 represents an often under- 
appreciated gain in sensitivity between a 21cm BAO experiment 
and a similar reionization experiment. At z = 9, X^Y ^ 
54:0(h~^Mpc)^ ST~^ Hz"-*^, meaning that the quoted EoR signal 
strength of ~ lOmK"^ is normalized over a much larger volume. 
The smaller volume scalar ai z = 1 means that over an order of 
magnitude less thermal noise is picked up per equivalent integra- 
tion. 
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in a considerably lower noise level, since our noise power 
spectrum scales as ft into equation [4j 

The issue of system temperature is another instance 
where a BAO experiment at z = 1 is fundamentally dif- 
ferent from an EoR experiment at 2; = 9. In the EoR 
case, Galactic synchrotron emission is the primary source 
of noise at 150 MHz, with a value of ~ 300K toward the 
galactic poles. However, synchrotron emission scales ap- 
proximately as 



<- sync 



300K 



(l50MHz) 



-2.5 



(11) 



At frequencies around 750 MHz, Tg, 



5K. This value 



is substantially below typical receiver temperatures of 
50K, so that receiver temperatures will dominate Tgyg. 
As described in ^ BAOBAB will have a system tem- 
perature of approximately 50K across its entire band. 

We can express the noise power spectrum of one dual- 
polarization baseline integrating on one Fourier mode by 
substituting these fiducial values in equation [4j 



A^(/c) ^ 8 X 10" 



0.1 hUpc 
180 days 



n 

0.045 sr 





'\u\- 




J-sys 




20 




50 K 



mK^ 



(12) 



where \u\ is the length of the baseline in wavelengths. 
The additional factors of sets the diameter of a uv- 
pixel and \u\ sets the length of time the baseline samples 
that same pixel. Our fiducial "long observation" is 180 
days; we set this value as a hard maximum for the num- 
ber of days BAOBAB can observe in one calendar year. 
This choice comes from the fact that observations will 
be compromised by foreground emission when either the 
Galactic plane or the sun is in view. 

In order to calculate the sensitivity of an entire inter- 
ferometric array, we must sum the contributions of all the 
Fourier modes sampled every baseline, paying careful at- 
tention to the number of times a unique Fourier mode is 
sampled by distinct redundant baselines. In our formal- 
ism, the effect of array configuration enters through the 
f / fo metric for the amount of redundancy in an array 
configuration, defined in P12a as: 



f_ 

fo 



En? 

i 



(13) 



where i labels individual uv pixels, and rii the number 
of one-second integration samples falling within pixel i. 
An ///o factor of 10^ is representative of the 32-element 
maximum redundancy arrays also described in PI 2a. 
The sensitivity boost comes from the fact that redun- 
dant sampling beats down thermal noise more quickly, 
since redundant measurements can be added in temper- 
ature, whereas distinct modes with the same |^| can only 
be added in the temperature-squared units of a power 



spectrum. The resultant sensitivity is: 



Af^{k) « 2 X 10" 



0.1 /iMpc"^ 



0.045 sr 



50 K 



■ 8 hrs " 


1 
2 


"180 days" 




"32" 




[lOVol 


_^per_day . 




^days 




_iV_ 




[ / I 



mK^. 

(14) 

The f / fo term is computed from uv sampling patterns 
including earth rotation aperture synthesis effects, and 
accounts for most all effects regarding array configura- 
tion. Hence, the additional factors of (7 2 and \u\ that 
appeared in equation [12] do not appear in[T4] However, 
we have separated the term containing the number of 
antennas, to give a sense of how sensitivity will roughly 
scale with array size. The factor of tper_day sets the to- 
tal integration time per day, which will likely be limited 
by the size of a low foreground emission region (i.e., a 
"cold patch"). We choose 8 hours as the maximum time 
that can spent observing cold patches per day, a value in- 
fluenced by existing all-sky maps and observations with 
PAPER. To calculate ///o, we split this 8-hour period 
into separate one-hour observations phased to a single 
pointing center, roughly corresponding to the width of 
the BAOBAB primary beam, after which a statistically 
independent patch of sky comes to dominate the data. 

3.2.2. Sample Variance 

In galaxy redshift surveys, the sample variance can be 
calculated relatively simply by counting Fourier modes 
over an effective survey volume. However, in the case of 
BAOBAB, not all modes are equal, since we have used 
redundant samples of certain modes to beat down ther- 
mal noise. In creating a 1- or 2-D power spectrum out 
of the full 3-D Fourier space, we must take a weighted 
combination of these modes, since the ratio of thermal 
noise to sample variance can vary between every k-mode 
measured. 

Using inverse- variance weighting to combine each mea- 
surement at a particular /c-mode, one finds that the opti- 
mal estimator of the power spectrum results in an error 
that can be calculated by combining the errors on each 
separate mode in inverse-quadrature: 



V 



AL(fc) 



^21 



{k)y 



(15) 



where 5A'^{k) is the resultant uncertainty on a given k- 
mode, is the thermal noise power spectrum calcu- 
lated in 



3.2.1 



A2i{k) is the cosmological 21cm power 
spectrum (which is also the sample variance error), and i 
is an index labeling the independent /c-modes over which 
we are summing (note that we never combine modes 
into bands with significantly varying k, which is why 
we can sum inverse absolute variance instead of inverse 
fractional variance). 

Since the sample variance is completely a function of 
the iii;-coverage of an array, it is best calculated numeri- 
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terms of 21cm brightness: 




- BAOBAB-32 Full Sensitivity 

- - BAOBAB-32 Sensitivity (Thermal Noise Only) 

BAOBAB-32 Sensitivity (No Foreground Cut) 

- BAOBAB-128 Full Sensitivity 

- . BAOBAB-128 Sensitivity (Thermal Moise Only) 

- ■ BAOBAB-128 Sensitivity (No Foreground Cut} 



k [h Mpc"^] 

Figure 8. Noise levels for two fiducial observations of one decli- 
nation range with BAOBAB at z = 0.89. The dashed thin (thick) 
line shows the sensitivity given thermal noise only for a 30 days ob- 
servation with a 35-element system (180 days with a 132-element 
system); solid lines show the effect of including sample variance. 
We also plot the effect of our foreground model of §3.3| the achiev- 
able sensitivity if no modes are excluded by foregrounds is shown 
by the dot-dashed lines. The other two curves do include the ef- 
fects of the foreground model. For comparison, we also plot the 
dimensionless 21cm power spectrum as the thick, gray line. At the 
larger scales, the 180 day observation with 132-tiles is completely 
dominated by sample variance. At the smaller scales, the analytic 
expression for thermal noise given in equation |14| accurately re- 
produces the sensitivity. The plotted thermal noise only curve is 
not a perfect power-law because longer baselines cannot probe the 
largest scale /c-modes. 

cahy. In Figure [8] we plot the ID thermal noise and sam- 
ple variance uncertainties for two maximum-redundancy 
configurations of BAOBAB — 35- and 132-tiles — shown 
in Figures [6] and [To) To calculate the sample variance, we 
use the 21cm brightness power spectrum from CAMB, 
also plotted for comparison. At the scale of the first 
acoustic peak, sample variance has clearly become the 
dominant source of error for a long integration with 132- 
elements; we discuss possible methods for improving this 
situation in §3.5| and §4.2[ At the smaller scales, how- 



ever, the analytic functions given in the previous section 
accurately reproduce the array sensitivity. Note that the 
thermal noise only curve in Figure [8] is not a perfect 
power law because not all baselines can probe the largest 
scales. We also plot the effect of our foreground model of 
§3.3] the achievable sensitivity if no modes are excluded 
by foregrounds is shown by the dot-dashed lines. The 
other two curves do include the effects of the foreground 
model. 

3.2.3. Shot Noise 

Measurements of the 21cm power spectrum will also 
be affected by the discrete nature of the neutral hydro- 
gen distribution at low redshift. Only overdensities self- 
shielded to the 1 Ry ionizing background contain neu- 
tral hydrogen, so we will be subject to the same galactic 
shot n oise as optical redshift surveys. Following Seo et al. 
( |2Q1Q| , we assume that galaxy positions and luminosities 
are distributed with probability proportional to l-\-b Sm^ 
where b is the bias and 5m is the mass overdensity. This 
allows us to treat shot noise as a scale-independent con- 
tribution to the power spectrum with P{k) = 1/n. In 



2 1 fc3 



(16) 



Using the result of |Seo et al.| (2010) (who use a fit to th e 
neutral hydrogen mass function of Zwaan et al. 2005), 

-3 



0.01 h^Mpc~ 
-1 



and A^^^, 



6.5 X 10-^mK^ at 
k = 0.1 hMpc~^ and z = 1. The number density of 
hydrogen- containing halos is substantially higher than 
for the bright galaxies used in optical/NIR surveys, mak- 
ing shot noise a substantially smaller contaminant; it will 
only begin to dominate the signal at /c > 2 hMpc~^ . Re- 
gardless of uncertainties in this calculation, shot noise is 
clearly a subdominant effect, and we neglect it for the 
remainder of this work. 

3.3. The Delay Spectrum Technique at z ^ 1 

Before combining the results of the last two sections, 
we must discuss the effect of foregrounds on power spec- 
trum measurements. The presence of foreground emis- 
sion orders of magnitude brighter than the cosmolog- 
ical 21cm signal has been one of the major impedi- 
ments for high-redshift 21cm tomography. PI 2b pre- 
sented a per-baseline delay-spectrum technique for isolat- 
ing foreground emission solely on the basis of its spectral 
smoothness. In this section we briefiy recapitulate the 
principles of the delay-spectrum technique, and present 
a simple approximation for the behavior of foregrounds 
in the 600-900 MHz band. 

The delay spectrum technique is a methodology for us- 
ing each baseline of an interferometer as a independent 
probe of the 21cm power spectrum. The most power- 
ful aspect of this approach is that the frequency depen- 
dence of a baseline's Fourier sampling pattern, typically 
regarded as a major complication for 21cm experiments, 
naturally gives rise to an isolation of foreground emission 
in Fourier space. The ability to remove foregrounds on a 
per-baseline basis allows multiple baselines to be tuned 
to target the same Fourier mode for greater sensitivity, 
as opposed to more traditional techniques that use over- 
lapping li'u- coverage at multiple frequencies to avoid the 
issue of frequency-dependent sampling. 

At the heart of the delay transform is a dual interpre- 
tation of the Fourier transform of interferometric visi- 
bilities along the frequency axis. On the one hand, for 
21cm experiments, frequency maps directly into redshift 
since the observed signal is a spectral line. Therefore, 
the Fourier transform along the frequency axis gives /cy , 
the Fourier wavemode along the line of sight. However, 
the frequency dependence of a baseline's length (as mea- 
sured in wavelengths), gives rise to the delay transform 
interpretation of the freque ncy Fo urier transform pre- 
sented in Parsons & Backer (2009). If performed over a 
wide enough bandwidth, this transform maps sources to 
Dirac delta functions in "delay space," corresponding to 
the geometric delay of signal arrival time between the two 
elements of the baseline. There is thus a maximum delay 
at which any signal coming from the sky can appear, set 
by the physical length of the baseline. Furthermore, each 
source delta-function will be convolved by a kernel rep- 
resenting the Fourier transform of that source's intrinsic 
spectrum (as well as any spectral features introduced by 
the instrument). Therefore, as long as the instrumen- 
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tal frequency structure is kept to a minimum, sources 
with smooth intrinsic spectra (such as foreground emis- 
sion) will have their emission confined within the region 
of delay space set by the maximum delays (the so-called 
"horizon-limit"). Sources with unsmooth emission, like 
the 21cm signal, will be convolved by a broad kernel, 
scattering "sidelobes" well beyond the horizon limit, and 
creating a window for detecting 21cm emission free of 
smooth-spectrum foregrounds. 

A full simulation quantifying the effects of the delay 
transform on foregrounds is beyond the scope of this 
work. Rather, we assign a baseline-dependent minimum 
below which we exclude modes as being wholly con- 
taminated by foregrounds. We chose the value of k 



based on the simu lations presented for PAPER in |Par- 
sons et al. (2012bl), which finds for the shortest baselines 
foregrounds contaminate modes with k\\ ^0.2 hMpc~ . 
At EoR frequencies of 150 MHz, this cutoff maps to 
delay-modes of 400 ns. Since BAOBAB baselines are 
physically shorter by a factor of 5, this reduces the max- 
imum delay-space contamination to 80 ns, which in turn 
maps back to /c = 0.1 hMpc~^ at 750 MHz, using the 
Y parameter from equation [S] There are two important 
factors which will further serve to reduce this number 
for BAOBAB. First, celestial foregrounds should have 
power law spectra with steeply decrease in intensity ver- 
sus frequency, and so will be fainter than at EoR fre- 
quencies. Although the signal has also fallen a similar 
amount, this reduced foreground structure will still make 
the delay transform even more effective at isolating fore- 
ground emission. Secondly, the narrower primary beam 
of BAOBAB will limit the delay modes from which there 
can be appreciable celestial emission, as sources near the 
horizon will be significantly attenuated. 

We perform a cursory calculation in which the delay- 
transform is applied to a simulated sky model. In these 
simulations, we assume the sky is entirely composed of 
point sources, where the source strength distribution fol- 
lows a power-law with a slope of -2.0, normalized to a 2 
Jy source per 10 steradians, a distribution derived em- 
pirically from PAPER data with extrapolation to the 
BAOBAB band. We also model the frequency spectrum 
of each source as a power-law with a normal distribu- 
tion of spectral indices centered on -1.0 and a standard 
deviation of 0.25. We refer to these simulations as "cur- 
sory," since they exclude instrumental effects such as RFI 
fiagging and frequency-dependent beam structure. In- 
stead, we use a single, frequency-independent Gaussian 
to model the primary beam of the BAOBAB tile; the 
potential effects of a more realistic beam model are dis- 
cussed in Q We find the delay transform confines fore- 
grounds to /c-modes below a value of k\\ = 0.045 /iMpc~^ 
for baselines of 16 wavelengths. The k\\ value for the 
maximum delay of a 16A baseline at 750 MHz (i.e., the 
horizon limit) is 0.028 /iMpc~^, which implies that the 
intrinsic spectral behavior of foreground emission corre- 
sponds to a kernel of width ~ 0.02 /iMpc~^. In this work, 
our foreground model is to exclude /c-modes smaller than 
the sum of the maximum realizable delay on a baseline 
(converted from seconds of light-travel time to hMpc~^ 
using equation [8| and this kernel. In effect, this model 
states that intrinsic spectral structure in foregrounds cor- 
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Figure 9. Predicted constraints on. z = 0.89 21cm power spec- 
trum from a 30 day observation with a 35-tile BAOBAB system. 
The net result is a ~ 5cr detection of the power spectrum. Results 
are comparable in the other 2 redshift bins. 

rupt 0.02 /iMpc~^ beyond a naive prediction based only 
on the physical length of the baseline. We explore the 
effects of modifying this model in §4.1| 

3.4. Detecting the HI Power Spectrum 

The first major science result from BAOBAB will be 
the detection of the 21cm power spectrum near 2: ~ 1. A 
high significance detection will be achievable with a short 
~ 1-month observation with a 35-tile system opera ting 
in the maximum redundancy shown in Figure |6p^| The 
predicted measurement for a 30-day observation (240 
hours) is shown in Figure |9] These observations assume 
a 100-MHz bandwidth centered on 750 MHz {z = 0.89). 
The net result is an b.6a detection of the 21cm power 
spectrum when our model for foreground emission from 
§3.31 is used to exclude contaminated modes. Results 
for bands centered on 650 {z = 1.18) and 850 MHz 
{z = 0.67) are similar, yielding 5.8<j and 5. Oct detec- 
tions, respectively. The lower significance detection at 
the lowest redshifts results from the scaling of the angu- 
lar diameter distance; at redshift z = 0.67, a 16A baseline 
corresponds to a k± wavemode of ~ 0.06 /iMpc~^, lim- 
iting the number of baselines that can probe the largest- 
scale /c- modes where thermal noise is lowest. Over the 
z = 0.5 — 1.5 range, both the 21cm signal and, the 
noise remain roughly constant in magnitude, the lat- 
ter because it is dominated by a frequency-independent 
front-end amplifier noise temperature. This trend does 
not continue indefinitely, however, as sky noise increases 
with increasing redshift, eventually dominating the sys- 
tem temperature. 

Measurements of this significance will allow for an ac- 
curate determination of /hi the combination of the 
cosmic neutral hydrogen fraction and the bias of neutral 
hydrogen containing regions, as a function of redshift. 
Breaking the degeneracy between these parameters will 
require additional information. Measuring redshift-space 
distortions can, in principle, separate the effects of the 

Correlator inputs have traditionally been in powers of 2; hence 
this array is our ~ 32-tile configuration. The sensitivity calcula- 
tions do assume that all 35 elements are correlated. 
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Figure 10. The maximum redundancy configuration of a 132- 
element BAOBAB system. The close-packed tiles are chosen to 
produce the shortest possible baselines. 

two terms. Constraints from a longer integration or a sys- 
tem with ~ 64 elements will further improve constraints 
on the neutral hydrogen power spectrum, and it will be 
possible for these systems to measure redshift-space dis- 
tortion effects. Measuring these effects accurately re- 
quires more careful systematic control, which may war- 
rant different configurations and observing strategies, so 
we postpone a n exploration of th is science to a future 
work (see, e.g., Masui et al.||201Q ). 



3.5. Detecting Baryon Acoustic Oscillations 

As shown in in Figure |8j a 132-tile BAOBAB array 
with the configuration shown in Figure [lO] has effec- 
tively reached the sample variance limit in 180 days (1440 
hours) of observing time. This observation yields a 3.3<j 
detection of the BAO features at z = 1.18, with effec- 
tive 2.1— and 2. 7 a non-detections at the z = 0.67 and 
z = 0.89 bands, respectively, where we have isolated the 
BAO features from the broad-band shape of the power 
spectrum by removing a model fit using the transfer func- 
tion from jEisenstein & Hui 



The effect of sample 
variance is most dominant at the lower redshifts, because 
the angular diameter distance scaling means that fewer 
samples of the BAO scale can be found in the same area 
of sky (this observation of 8-independent fields with a 
0.045 steradian primary beam corresponds to an effec- 
tive survey area of ~ 1200 square degrees). While longer 
observations with the same array configuration can im- 
prove these constraints by reducing thermal noise on the 
smaller scale modes, a better approach will be to observe 
additional independent fields. 

For our fiducial BAOBAB observation, we use an ar- 
ray which observes 24 independent fields (i.e., three in- 
dependent declination observations for 8 hours per day), 
yielding an effective survey area of ~ 3600 square de- 
grees. We discuss the motivation for thi s pa rticular ap- 
proach to increasing survey volume in ^4.2 There are 
several equivalent ways an experiment can probe this ad- 
ditional area. One approach would be to conduct three 
1-year surveys, with the dipoles pointed towards a dif- 
ferent declination; this could be achieved by physically 
placing the array on an platform inclined by ~ 15°, or 
potentially by adding a steering component to the tile- 
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Figure 11. Predicted constraints on the z = 0.89 BAO features 
from a 180 day observation of three-declination fields with a 132- 
element BAOBAB system. The net result is a 4.7cr detection of 
the BAO features. The sensitivity is comparable in the other 2 
redshift bins. 

beamformers. If the beamformers are designed to allow 
multiple beams, one could in principle achieve similar 
sensitivities with only one year of observing, although at 
the expense of additional degrees of complexity in the 
system. If funding permits, the simplest approach might 
be to build three 132-tile BAOBAB arrays, each tilted 
towards a different declination; this would also yield the 
subsequently predicted sensitivities in only one season 
of observing. Note that for BAO science these indepen- 
dent configurations are potentially more desirable than 
an array with a larger number of fully correlated tiles. 
Since we are using a close-packed configuration, the ad- 
dition of more tiles can only yield new modes at corre- 
sponding larger k± where the amount of BAO informa- 
tion is significantly diminished. It may be possible that 
more information could be recovered from these larger 
k modes using a reconstruc tion method (e.g. Eisenstein] 



et al. 2007, Padmanabhan et a! 



Noh et al I I 



but we do not explore this option m this work. It there 
is significant BAO information beyond k ~ 0.2, then a 
~ 256-tile array could possibly yield tighter constraints 
than two 132 -tile arrays observing for the same amount 
of time. In §4.2[ we discuss other ways to increase the 
survey footprint and further reduce the effect of sample 
variance. 

We plot the expected constraints for the three- 
declination range fiducial BAOBAB observations on the 
z = 0.89 BAO features in Figure [TT] where the broad 
band shape of the power spectru m has been removed us - 



ing the transfer function fit of Eisenstein & Hu ( 1998 ) 
which neglects BAO. The measurements from a 180 
day integration at each declination range with this ar- 
ray amount to an 4.7 a detection of these features. Re- 
sults are similar for the other redshift bins, with expected 
3.6(7 and b.7a detections at redshifts of 0.67 and 1.18, re- 
spectively. While it is clear from Figure |8] that sample 
variance dominates the errors on the largest scale modes 
after 180 days of observing one declination range, we find 
the additional sensitivity towards the higher BAO peaks 
with this observing duration yields better constraints on 
the signal than e.g., observing twice as many declination 
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ranges for 90 days. 

With a significant BAO detection, BAOBAB can also 
begin to place constraints on cosmological parameters. 
To quantify the effect of such measurements, we use the 
Fisher matrix formalism of the Joint Dark Energy Mis- 
sion ( JDE M) Figure of Merit Science Working Group 
(FoMSWG; |Albrecht et al.|2QQ9D , defining our Fisher ma- 
trix as: 



(17) 



where / is some observable measured at some b values, 
cr^ is the variance in a measured value of /, are cos- 
mological parameters, and we sum over all measured 
values. 

We propagate our power spectrum measurements into 
constraints on the Hubble parameter H{z) and angu- 
lar diameter distance Da using this formalism, where 
/ = A^(/c), our measured power spectrum, b = the 
set of /c-modes we measured, and = [H{z)^ Da]- The 
derivatives with respect to Da and H{z) are straight- 
forward to calculate, as they affect our measurements 
through the X and Y parameters defined in equations 
[s] and [s] In effect, changing H{z) or Da changes the 
T- modes sampled by BAOBAB. For our calculations, 
we exclude fe-mo des d eemed contaminated by our fore- 
ground model of §3.3[ We also model the nonlinear de- 
generation of high er fe-modes using the elli ptical Gaus- 
sian formula from Seo & Eisenstein (2007). To isolate 
the constraints provided by the BAO features from the 
broad-band shape of the power spectrum, we again re- 
move a model fit u sing the transfer function from Eisen- 



stein & Hu (1998). We split our data into three redshift 
bins centered at z = 0.67,0.89 and 1.18. Although our 
frequency coverage is continuous between z = 0.58 and 
1.37, we find that there is minimal penalty for using only 
three bins in a Fisher matrix study. The result of this cal- 
culation is that our fiducial 3-declination, 180-day inte- 
gration yields measurements of H{z) with an error rang- 
ing from 9% to 4.5% across our three redshift bins, from 
low to high redshift, and measurements of Da with er- 
rors effectively constant at 17% over the same range (note 
that errors > 5% should be understood in the usual for- 
mal Fisher matrix sense - these measurements would not 
correspond to significant detections, taken in isolation). 
The particularly poor constraints on Da come from the 
loss of modes due to foreground emission; w e fur ther ex- 
plore the effects of our foreground model in ^4.1 Taking 
correlations between H{z) and Da into account, these 
measurements amount to 2.5% to 4.5% errors on a "dila- 
tion factor," which scales Da and H~^{z) in proportion. 
Expressed as a single constraint on the z ^ 1 distance 
scale, these measurements correspond to a dilation factor 
error of 1.8%. The exact uncertainties are given in Table 

m 

It is also straightforward to propagate these errors on 
H{z) and Da into errors on underlying cosmological pa- 
rameters through the Fisher matrix formalism of equa- 
tion [TtI In this case, / is now H{z) or Da^ and b cor- 
responds to the redshift bin it was measured in. H{z) is 



given by: 

H\z) = Hi + zf + + zf^ 

^D£;exp 3 / — ^[l+w(zO] 
Da is given by: 



cdz 



(18) 



(19) 



The parameters of interest, p% are now the underlying 
cosmological parameters. We use the parameterization 
of the JDEM FoMSWG, which include Wm^Wb^Wk and 
wde {u)x = ^xh^^ where m^b^k^ and DE correspond 
to the matter, baryon, curvature and present day dark 
energy density, respectively) although we use the simpler 
2-component form for the dark energy equation of state: 



w{a) = + (1 - tt)^a- 



(20) 



Follo wing the convention o f Dark Energy Task Force re- 
port ( Albrecht et al.|2006 ), we marginalize over all other 
parameters alter conibining our Fisher matrices with con- 
straints from other experiments, to create 2x2 matrix 
representing constraints on wq and Wa- As a Figure of 
Merit we use FoM = where J^' is our original 

Fisher matrix J^, marginalized to a 2D (k;o, '^a)-space; 
this FoM is proportional to the inverse of the error el- 
lipse area in the — Wa plane. 

As a baseline for current dark energy constraints, we 
use the JDEM FoMSWG predictions for the Planck satel- 
lite, combined with constraints from the BOSS-LRG sur- 
vey listed in Table [2] and a 5% error on Hq. This com- 
bination of experiments yields an FoM of 8.7. Including 
the measurements from our fiducial BAOBAB observa- 
tion increases this FoM value to 16.6. This FoM shows 
that the high-redshift constraints from BAOBAB can sig- 
nificantly improve our current measurements of the dark 
energy equation of state. These constraints are plotted 
as 1 and 2cr error ellipses in Figure [121 If we include 
the Lya forest survey of BOSS, our baseline constraint 
FoM becomes 17.8 which is improved to 23.0 with the 
inclusion of BAOBAB data. Even when the BOSS-Lya 
forest constraints between redshifts 2 and 3 are added, 
BAOBAB still provides valuable information, serving to 
increase the FoM by - 25%. 

4. DISCUSSION 

We break this discussion into two parts. In §4.1[ we 
consider two components of our analysis which may be 
overly simplistic: our model of the BAOBAB primary 
beam as a frequency-independent Gaussian, and our fore- 
ground emission model. We discuss the effects any short- 
comings in these models could have on our conclusions. 
In §4.2[ we consider the fact that sample variance is the 
dominant source of uncertainty in our measurements, 
and present future approaches that could improve the 
dark energy constraints possible with the BAOBAB in- 
strument. 

4.1. Potential Shortcomings in the Analysis 

A effect that could modify the predictions presented in 
this work is the use of an overly simplistic model for the 
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Table 2 

Percent errors on the distance scale from BAO measurements, for a three dechnation BAOBAB survey. The correlation is the correlation 
coefficient between the H(z) and measurements. R is the "dilation factor," a single estimate of the distance scale which scales 
and H~^{z) in proportion. BOSS-LRG constraints come from ^ Schlegel et al. (200 9| and BOSS-Lya Forest constraints come from the 

method of jMcDonald Eisenstem, ( ,:^(JU7^ ) and |Dawson et al.] ( |2012| >. 



Survey 


Redshift 


H{z) Error 


Da Error 


Correlation 


R Error 


BAOBAB 


0.67 


8.9% 


17.1% 


0.71 


4.4% 


BAOBAB 


0.89 


6.1% 


16.4% 


0.72 


3.3% 


BAOBAB 


1.18 


4.5% 


17.5% 


0.73 


2.6% 


BOSS-LRG 


0.35 


1.8% 


1.0% 


0.41 


0.7% 


BOSS-LRG 


0.6 


1.7% 


1.0% 


0.41 


0.7% 


BOSS-Lya 


2.5 


3.1% 


7.4% 


0.58 


2.0% 



— BOSS-LRG + Planck + BAOBAB 

- BOSS-LRG + Planck 




-1.(1 



Figure 12. 1 and 2cr error ellipses in the WQ—Wa plane for various 
surveys. The dotted line shows constraints from Planck, BOSS- 
LRG and a 5% error on Hq] the solid line shows the effect of in- 
cluding a 1440-hour integration on three independent declination 
fields with a 132-tile BAOBAB array. These error ellipses corre- 
spond to Figures of Merit of 8.7 and 16.6, respectively. 



BAOBAB beam: a frequency independent Gaussian. Al- 
though this model is sufficient to calculate sensitivities, 
the effects of a more realistic beam model on the delay 
spectrum foreground removal technique will necessitate 
further investigations. The principle cause for concern is 
from the existence of frequency- dependent grating-lobe 
structure associated with the tile. These sidelobes can 
introduce apparent frequency structure into otherwise 
spectrally smooth foreground emission. If this structure 
represents a significant increase in the size of the delay- 
space convolving kernel, it will move the foreground con- 
taminated region to higher k values. If the size of the 
effect is large enough to push foreground emission on the 
shortest baselines beyond the first BAO peak, the pre- 
dicted cosmological constraints could be reduced. 

Note, however, that the frequency-independent Gaus- 
sian beam is not as bad an assumption as it might first 
appear. If we were trying to image the sky, ignoring the 
grating lobes associated with the tile would be unjus- 
tified. With the delay spectrum approach, though, the 
issue is not the existence of the sidelobes, but their fre- 
quency dependence. If the frequency Fourier transform 
of the beam pattern is particularly broad — correspond- 
ing to rapid evolution of the beam pattern with frequency 
— then foreground emission will have a similarly broad 
footprint in delay space, compromising the 21cm signal. 
However, as argued in PI 2b, it is difficult for an element 
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Figure 13. The sampling pattern of BAOBAB-132 at ^ = 0.89 in 
the (fc_L-fc|| )-plane. Modes are marked as either uncontaminated or 
contaminated by foregrounds using our model of §3.3| Foregrounds 
limit jji to, e.g., > 0.65 for k ~ 0.1, leading to the relatively poor 
Da measurement and high correlation. 

only several wavelengths across to possess such frequency 
structure. As stated in §3.3[ we will postpone a detailed 
investigation of these effects to a future paper, and focus 
here on array sensitivities, for which the Gaussian model 
is sufficient. 

As noted in §3.5[ our constraints on Da are signifi- 
cantly limited by foreground emission. We illustrate this 
effect in Figure [Tsj which shows the sampling pattern of 
BAOBAB-132 at 750 MHz in the (A:^-A:||)-plane, high- 
lighting those modes discarded as contaminated by our 
foreground model of ^3.3 Foreground contamination ef- 
fectively excludes modes where < /cy, i.e., transverse 
modes. This has the effect of significantly degrading the 
achievable constraints on Da- 

While our foreground model presented in §3.3| is em- 
pirically motivated, accurate predictions for foreground 
emission will need to come from early BAOBAB arrays or 
other 21cm experiments. We briefly explore the effect of 
changing the size of the additive component of our fore- 
ground model arising from the spectral smoothness (or 
lack thereof) the emission. For our fiducial model, this 
term has a magnitude of k\\ = 0.02 /iMpc~^; as test cases, 
we analyze the constraints obtainable with BAOBAB if 
this term is changed by A/cy = ± 0.01 /iMpc~^. This 
has the effect of moving the cutoff between contaminated 
and uncontaminated modes in Figure [13] up or down by 
0.01 /iMpc"\ 

The effect of increasing this term (i.e., if foreground 
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emission is not as spectrally smooth as predicted) is to 
degrade the significance of our BAO detections in each 
of the redshift bins by ~ 0.4cr, our H{z), Da^ and R 
constraints by ~ 0.5, ~ 5, and ~ 0.3 percentage points, 
respectively. Reducing the foreground emission footprint 
in /c-space (i.e., if foregrounds are spectrally smoother 
than predicted) has similar effects with the opposite sign: 
the significance of our BAO detections in each band are 
increased by ~ 0.3cr, our H{z)^ Da, and R constraints are 
improved by ~ 0.5, ~ 5, and ~ 0.3 percentage points, re- 
spectively. From this analysis, it is clear that foreground 
emission can significantly alter the achievable constraints 
on Da, but ultimately, the success of BAOBAB will not 
be determined by the details of foreground spectral prop- 
erties. 

4.2. Improving Our Constraints 



As discussed in §3.5[ sample variance dominates the 
uncertainties in power spectrum measurements from 
BAOBAB. However, measuring a new set of indepen- 
dent modes is not a trivial undertaking. By using a 
close-packed array, BAOBAB completely samples the uv- 
plane out to some maximum baseline length. Therefore 
a simple array rotation or reconfiguration will not yield 
new samples. Furthermore, in our highest redshift bin 
of z = 1.18, the longest baseline in the 132-tile array 
probes a transverse /c-mode of k± = 0.15 hMpc~^ (and 
at lower redshift, this longest baseline corresponds to an 
even larger value of k±). At this scale and smaller, most 
BAO information is being lost to non-linear damping 
effects. Therefore, the returns from increasing the ar- 
ray size beyond ~ 128 tiles are limited, since effectively 
no new modes with significant BAO information will be 
probed. 

As an upper-limit to the constraints obtainable with 
a single-declination, ~ 1200 square degree, 132-tile 
BAOBAB observation (as opposed to our fiducial ob- 
servation targeting three declination fields) we compute 
the results of a completely sample variance limited sur- 
vey, i.e., one where thermal noise uncertainties have been 
set to (although modes are still excluded using our 
foreground emission model). This sample variance lim- 
ited 132-tile BAOBAB observation yields distance scale 
uncertainty of ~ 1.5%, averaged over our whole band 
(compared with 1.8% for our fiducial 1-year observa- 
tion of 3 independent declination fields). For compari- 
son, a 10-year (14,400 hour), one declination observation 
with the same array yields a distance scale uncertainty 
of 1.9%. Therefore, while better constraints can come 
from a longer observation, obtaining measurements of 
new modes to beat down sample variance is clearly the 
optimal way to proceed. There are two ways forward to 
achieve this goal: map a different volume of the universe 
(as with our fiducial experiment) or recover foreground 
corrupted modes. We consider each of these approaches 
in turn. 

Since BAOBAB is a zenith-pointing, drift-scanning 
telescope, to map a new area of the sky, it will either 
need to be relocated to a different latitude or tilted to 
point towards a different patch of sky. Either option is 
potentially feasible, as even a 132-tile BAOBAB experi- 
ment spans less than 20 meters. With a primary beam 
full- width half-max of ~ 15°, there are ~ 10 independent 
pointings in declination that BAOBAB can target. Since 



our fiducial observation already targets three declination 
fields, mapping every declination could in principle yield 
up to ~ 70% reductions in the error bars over the results 
presented. 

The other way to potentially measure new modes and 
beat down sample variance is to recover samples we have 
considered corrupted by foregrounds. There are two ways 
foregrounds compromise BAOBAB observations. The 
first is the limited observing time per-day, set by Galactic 
emission, which we have treated as irreparably corrupt- 
ing all samples, even those in principle recoverable with 
the delay transform. If it is possible to observe all 24 
hours of right ascension, as opposed to the 8 considered 
here, the constraints from a single observing season will 
increase by a factor of a/3. While it is unlikely that all 
24 hours of right ascension will be workable, our fiducial 
value of 8 hours per day, motivated by observations with 
present EoR experiments, may well be conservative, since 
Galactic synchrotron emission has significantly fallen in 
brightness compared to EoR frequencies. 

Even when observing a "cold patch," foregrounds cor- 
rupt large scale /c-modes with a footpri nt m oving to 
smaller scales as baseline length increases (^3.3). If these 
modes could be retrieved, they could signiticantly in- 
crease the volume of Fourier space that BAOBAB can 
probe. As an upper-limit, we calculate the obtainable 
power spectrum constraints ignoring all foreground con- 
tamination. The result is that a 1-year (1440- hour) ob- 
servation in each of three independent declination fields 
yields a distance scale uncertainty of 1.4% combined over 
the entire redshift range, an increase of ~ 25% over the 
same observation including foreground emission. In par- 
ticular, we note that a foreground- free observation yields 
errors of - 5% on Da at redshifts of 0.67, 0.89 and 1.18, 
respectively — a factor of > 3 improvement over the 
predictions for an observation including the effects of 
foreground emission. While an analysis of foreground 
removal techniques is beyond the scope of this work, this 
result is suggests that foreground removal may be the 
way to improve constraints on Da- 

As an order of magnitude estimate, we can consider 
whether a foreground removal or subtraction scheme 
might be more effective in the BAO band than at EoR 
frequencies. At /c - 0.1 /iMpc"\ the 750-MHz BAO 
21cm power spectrum reaches ~ 3 x 10~^ mK^, com- 
pared with a peak brightness at EoR frequencies of 150 
MHz reaching ~ 10 mK^. The steep spectrum Galactic 
synchrotron emission has a spectral index of -2.5, and so 
win fan by a factor of {5-^-^)^ = 3.2 x lO""^ in units of 
temperature squared. (Extragalactic point sources are 
less steep spectrum, and so will not fall off in bright- 
ness as steeply. Therefore, this estimate can be con- 
sidered a lower limit on the foreground-to-signal ratio). 
Roughly speaking, then, the foreground-to-signal ratio 
is unchanged compared with EoR experiments, suggest- 
ing that a foreground isolation scheme like the delay- 
spectrum technique is still likely the most viable ap- 
proach for first-generation experiments limited in collect- 
ing area. 

5. CONCLUSIONS 

In this work we have presented a concept for a new 
experiment using the redshifted 21cm line of neutral hy- 
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drogen to probe cosmology at 2; ~ 1. The BAO Broad- 
band and Broad-beam Array (BAOBAB) will incorpo- 
rate both the hardware and analysis infrastructure de- 
veloped for 21cm experiments at higher redshifts. 

The hardware design will borrow heavily from the Pre- 
cision Array for Probing the Epoch of Reionization (PA- 
PER) and the Murchsion Widefield Array (MWA), us- 
ing a scaled version of the PAPER dipoles as a feed el- 
ement, tiling dipoles as done by the MWA, and modi- 
fying the CASPER FPGA/GPU PAPER correlator to 
perform full dual-polarization cross-correlations of all el- 
ements. Significant improvements to the system tem- 
perature will be brought about through state-of-the-art 
uncooled, low-noise amplifiers. Relative to PAPER, the 
collecting area will be substantially increased through 
the use of tiles of 4 dipoles combined through a beam- 
former, as demonstrated by MWA efforts. Although we 
have largely avoided specific cost-estimates, it is fair to 
say that this infrastructure is obtainable at a fraction of 
the cost of the ~ 100-million dollar ground-based optical 
redshift surveys. 

On the analysis side, BAOBAB will use the max- 
imum redundancy configurations and delay spectrum 
foreground removal te chniques prese nted in Par sons et al.| 
|2012a) and Parsons et al. (2012b) to enhance sensitiv- 
ity to Fourier modes along the line-of-sight. Motivated 
by the science of Baryon Acoustic Oscillations, BAOBAB 
will utilize extremely close-packed arrays to maximize the 
number of short baselines. The sensitivity calculations 
presented here show that BAOBAB will achieve several 
milestone measurements over our anticipated staged de- 
ployment process. A ~ 32-element BAOBAB system will 
yield high significance detections of the HI power spec- 
trum, and constrain the evolution of the cosmic neutral 
hydrogen fraction from z = 0.5 to 1.5 as well as the 
bias of DLAs. Over the same wide redshift range, a 
~ 128-element system will allow for a first detection of 
the BAO features in the power spectrum, and yield er- 
rors on the distance scale R at the several percent level. 
When combined with our current constraints on dark 
energy, including those forthcoming from the BOSS and 
Planck experiments, BAOBAB's measurements result in 
substantial increase in the Dark Energy Task Force Fig- 
ure of Merit, representing constraints on the nature and 
time evolution of dark energy over a wide range of cosmic 
history. 
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